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In the present paper we study chemical potential of the generalized Hubbard model with corre- 
lated hopping. The peculiarity of the model in comparison with similar generalized Hubbard models 
is the concentration dependence of hopping integrals. Chemical potential as a function of the model 
energy parameters, electron concentration and temperature is found. It is shown that correlated 
hopping and temperature changes essentially the chemical potential location; these dependencies 
C ' differ strongly at different values of the electron concenntration. 
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I. INTRODUCTION 



One of the simplest models, which is used for a description of strongly correlated electron systems (in particular, 
transition metal compounds, polymers, fuUerenes Ceo, high-temperature superconductors, heavy fermion substances) 
is the Hubbard model (for reviews see papers Q). Hubbard first showed that taking into account the intra- 
atomic Coulomb repulsion U of two electrons essentially modifies the band energy spectrum and plays the main 
role in the formation of electrical and magnetic properties of the material with narrow energy bands. The ideas of 
d , "lower" and "upper" Hubbard bands, introduced in this paper, proved to be useful both for our understanding of 
physics of correlation effects in narrow energy bands and for the interpretation of experimental data for narrow band 
I materials, in particular high- Tc superconductors. It should be noted that progress in experimental techniques, such as 
' O photoemission and x-ray absorption spectroscopy, has recently confirmed and enabled direct measurements of these 
^ ] theoretically predicted Hubbard bands jsj- 0. 

However, theoretical analyses, on the one hand, and available experimental data, which argue an electron-hole 
asymmetry, on the other hand, point out the necessity of the Hubbard model generalization by taking into account 
correlated hopping (electron- hole asymmetry is a peculiar property of the Hubbard model). Really, in some compounds 
(e.g. see the estimation in Refs. [Q) the matrix elements of electron-electron interactions describing inter-site 
hoppings of electrons (correlated hopping) arc of the same order that the hopping integral or on-site Coulomb repulsion; 
also, the results of the first principles calculations performed by Hirsch show that hopping integrals of electrons 
' between neighbouring ions can depend strongly on the instantaneous charge occupation of the two atoms involved 
\^ : in the hopping process (in contrast to the Hubbard model where inter-site hoppings of electrons do not correlate 
by the occupation of these sites). On the other hand, the experimental observations show electron- hole asymmetry 
, of properties of strongly correlated electron systems. In particular, transition metal oxides exhibit the electron-hole 
■ asymmetry of conductivity |]l6|, transition 3c?-metals have the electron- hole asymmetry of cohesive energy p^ , p^ . 
In this connexion, the experimentally observed |l9| , ^0t] electron-hole asymmetry of superconducting properties of 
high-temperature superconductors should be also noted. 
I Consequently, in recent years the generalized Hubbard model with correlated hopping has been used widely to 
. describe strongly correlated electron systems p^ ; the electron-hole asymmetry is a property of such a generalized 
(~| Hubbard model as a result of the dependence of the hopping integral on the occupation of the sites involved in the 
O hopping process. 

O Generalized Hubbard model has been used for an investigation of metallic ferromagnetism in narrow energy 

^ [ bands P,|2l|]- l2^ . In particular, a generalization of Nagaoka's theorem has been proved in paper and it has 
been shown p,p3[ that in strong coupling regime the model strongly favours ferromagnetism for electron concentra- 
tion n > 1 versus electron concentration n < 1, as observed in transition metal compounds. 

Electrical properties, metal-insulator transition and high-temperature superconductivity have been studied in the 
" " ' generalized Hubbard model in papers |jlJ,Q- The result of these papers is the generalized Hubbard model has 
much better physics than the Hubbard model, firstly, and usage of the electron-hole asymmetry concept allows one to 
interpret the peculiarities of physical properties of strongly correlated electron systems, which are not explained by 
the Hubbard model, secondly. In particular, the experimentally observed electron-hole asymmetries of metal oxides 
conductivity, of cohesive energy of transition 3d-metals and of superconducting properties of high-temperature super- 
conductors, mentioned above, have been explained within the generalized Hubbard model with correlated hopping in 
papers IIItIII, II, II respectively. 

An important puzzle remaining beyond consideration in the cited papers is the question about location of the 
model chemical potential and its dependence on electron concentration and temperature. Theoretical calculations of 
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chemical potential dependences can be compared with the experimental data obtained by photoemission spectroscopy 
methods. Also these dependences of chemical potential could give indication for the existence of phase separation in 
system upon doping. 

Even in the framework of the Hubbard model this problem is still open, in particular the question about the jump 
in chemical potential of the Hubbard model at half-filling in going from electron to hole doping is under discussion 
in recent literature (see, for instance, review [^). On the one hand, in majority of papers it has been found that 
the jump in chemical potential at half-filling is equal to width of the energy gap. On the other hand, some authors 
suggest that the discontinuity in chemical potential is strictly smaller than the energy gap as a result of inducing 
states inside the Mott- Hubbard gap (so called "midgap states" ) upon doping the Mott insulators. The first scenario 
is confirmed by the exact solution |371| of the Hubbard model in one dimension. Also numerical studies of the two- 
dimensional Hubbard model performed in Rcf. fsM] by means of the Monte Carlo methods show that when the electron 
concentration n — 1 — 6 {6 is infinitesimal doping) , the chemical potential is located at the top of the lower Hubbard 
band; when the electron concentration n = 1 + d, the chemical potential jumps to the bottom of the upper Hubbard 
band. Correspondingly, upon electron doping the chemical potential jumps from the top of the lower Hubbard band 
to the bottom of the upper Hubbard band according to the first scenario. The results of papers |39y40f| obtained by 
means of the Kotliar-Ruckenstein slave bosons p| for the three dimensional Hubbard model and the quantum Monte 
Carlo method in infinite dimensions suggest that the jump in chemical potential at half-filling is equal to width 
of the energy gap. 

Using the projective self-consistent technique in infinite dimensions, the authors of papers jij] found that, for 
any value of U larger than Uc, doping does induce states inside the Mott-Hubbard gap; the splitting of these states 
from the edge of the Hubbard band is a fraction of the bare kinetic energy, and remains finite even in the limit where 
U is infinite. Therefore the jump in chemical potential for infinitesimal doping is strictly less than the energy gap, 
so the second scenario described above holds. These conclusions were also supported by an extention of the 

iterated perturbative theory approximation |^] for the Hubbard model away from half-filling. Midgap states have also 
been predicted |^ to occur upon doping in the infinite dimensional Hubbard model using an analytic variant of the 
Lanczos continued fraction method |48|1 . 

In consequence of the absence of a rigorous result on the described scenarios adequacy the subject is still open. 
Moreover, even in infinite dimensions this problem is under discussion: different approximations, for example the 
quantum Monte Carlo method and the iterated perturbative theory, lead to the contrary results. Another question 
is how correlated hopping does infiuence on chemical potential. Therefore, the investigation of chemical potential 
behaviour in the doped Hubbard model as well as in the generalized Hubbard model with correlated hopping is a 
relevant and important task, which the present paper is devoted. 

The paper has the following structure. In Section 2 we formulate Hamiltonian of the generalized Hubbard model 
with correlated hopping, the quasiparticle energy spectrum of the model is derived by means of a generalized mean- 
field approximation (an analogue of the projection operation) in the Green function technique. The general properties 
of the spectrum, and the approximation used by us are discussed. In Sections 3-4 the dependences of chemical 
potential on the model energy parameters, correlated hopping, temperature and electron concentration are studied. 
The obtained results are compared with the experimental data on Fermi energy behaviour of real strongly correlated 
electron systems. Finally, Section 5 is devoted to the conclusions from the obtained results. 



We start from the natural generalization of the Hubbard model including the matrix elements of electron-electron 
interaction 



which correlate the electron hopping from one lattice site to another {(p{r — R^) is the Wannier function). The 
Hamiltonian of such a generalized Hubbard model with correlated hopping reads as [p4|,p6[ 



II. HAMILTONIAN AND QUASIPARTICLE ENERGY SPECTRUM OF THE MODEL 




(2.1) 
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H[=t'{n)} (Tj,XfXf + h.c 



ija 



(2.5) 



where /i is the chemical potential, the primes at the sums sig nify that j, = is the Hubbard operator [Q, 
X^ — X^^Xl^ is the operator of the number of |/c)-states on i-site; |0) denotes the state of site, which is not occupied 
by an electron (hole), \(t) denotes the state of singly occupied (by an electron with spin a) i-site, |2) denotes the state 
of doubly occupied (by two electrons with the opposite spins) i-site (doublon), t^i = — 1, y^j = 1. 



t{7i) = to + '^^^ J{ikjk) ~ to + nTi, 



i{n) = t{n) + 2T2, t' (n) = t{n) + T2, T2 = J{iiji) = J{iiij), 



(2.6) 
(2.7) 



where to is the uncorrelated hopping integral (the matrix element of electron-ion interaction), n ~ Ne/N is the 
electron concentration {N^ is the number of electrons within conduction band, N is the number of lattice sites). 

Hq describes the atomic limit of narrow-band models. Hi describes the translational hopping of holes and doublons; 
within the present model the hopping integrals of holes t{n) and doublons t{n) are different (in contrast with narrow- 
band models of the Hubbard type), this leads to the electron-hole asymmetry mentioned above. iJ{ describes the 
processes of paired creation and destruction of holes and doublons. 

In the model described by Hamiltonian ( ^.2[ ) an electron hopping from one site to another is correlated both by 
the occupation of the sites involved in the hopping process (with the hopping integral T2) and the occupation of 
the nearest-neighbour sites (wit h th e hopping integral Ti) which we t ook into account by means of the Hartree-Fock 
approximation (see expression ( [2.(j )). The peculiarity of model (2.2) is taking into consideration of the correlated 
hopping Ti which cause the concentration dependence of the hopping integrals in contrast with similar generalized 
Hubbard models. 

For half-filling case we introduce the following notations: tjn)^ t = to +Ti, i{n) = i — t + 2T2, t'{n) = t' = t + T2. 
In this case for t' = some exact results have been found |2|,|2^,^,Q. In a simple cubic lattice with coordination 
number z metal-insulator transition occurs at 



Uc = 2z\to\ 



(2.8) 



li U > Uc the ground state of system is a paramagnetic Mott-Hubbard insulator with the concentration of doubly 
occupied sites d = 0, the ground state energy is equal to zero. 

The single-particle Green function in terms of the Hubbard operators reads as 



= {{x'^^ix;^)) + vA{x"/\x'/)) + vA{xp\x;P)) + {{x;'\x^/)). 

The Green function {{Xp'^\X^P)) is given by the equation 

Sn 



(2.9) 



{E + ^,){{x^>^\x;P)) 



^^.{x; + x^,) + {{[x'/,Hi] \x;P)) 



(2.10) 



with \A, B ] = AB — BA. Using a variant 1 24 2^ of the generalized mean- field approximation |^0[ we suppose in 
Eq. (|2.10D that 



(2.11) 



where e (pj) a nd ei(pj) are the non-operator expressions. The procedure of ejpj) and ei{pj) calculation is described 
in Ref. l2J,Eq| (here there is a partial equivalence with the slave boson method pl[). 

Thus we obtain the closed system of equations for the Green functions {{X2^\Xy)) and {{Xp^\XpP)). An analogous 
procedure is realized also in the equations for the other Green functions ( ^.Q] ). 

In this way, we find single-particle Green function (2.9) and quasiparticle energy spectrum. For paramagnetic case 
in k-representation the spectrum is [p7|: 



/, ^ U e(k)-f-?(k) 1 , 

0(k) = v/[e(k)-g(k)-C/]2+4ei(k)e2(k). 



(2.12) 
(2.13) 
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The Fourier components of the quantities defining quasiparticle energy spectrum ( ^.12 ) are given by the formulae: 

e(k) = atkH, ?(k) = aikin), ei(k) = ai4(n), e2(k) = a24(n), (2.14) 



2d 



a = n — 2d ' 



2(l-rf)2 2d{d - n + I) i{n) 
2 - n 2~n t{n}' 

2d^ 2d(d -n+1) t(n) 
n n t{n) ' 

2d _ 2(1 - d) 

2-n ' 



ai = n — \ , q;2 = — 1 — ?i 

with d beeing the concentration of doubly occupied sites. 
For n < 1 and [/ ^ oo we obtain 



(2.15) 
(2.16) 
(2.17) 



i^i(k) = -M 



2-n 



- n tk(n) 



(the lower Hubbard band); if n > 1 and only the upper Hubbard band is important, we obtain 

E2(k) = -fi + U + (^ -2 + nA tk(n). 



(2.18) 



(2.19) 



For the Hubbard mod el (t{ n) — i{n) — t'{n) ~ to) -Ei(k) for n ^ 0, and -E2(k) for n ^ 2 get the band form. 
Moreover, expressions (2.12) describe the exact atom ic limit for to — and the band situation for U = 0. 
At half- filling qausiparticle energy spectrum ( 2.12| ) is writen in the form 26 1: 



i^k = V [B{tk - <k) -U] + [Mt'J^, B = l-2d + 4^2, 
where the doublon concentration d is found from the equation 



exp 



1 exp 



with 



. ^ 1 _ B{i^ - fk) _ _U_ 
^ 2 2Fu 2Fk' 



Bk = - 



1 , S(fk-ik) 



2Fk 



U 
2K' 



(2.20) 
(2.21) 

(2.22) 
(2.23) 



6 = ksT, ks is th e Bol tzmann's constant, T is the temperature . Fo r the special case t' = of the model quasiparticle 
energy spectrum ( 2.20| ) reproduces (see Ref. |^) exact result (2^ 



As a result of the absence of a natural small expansion parameter at intermediate to strong Coulomb interactions 
in the Hubbard model and its generalizations (this situation is under consideration in the present paper), one does 
not know a rigorously justified nonperturbative approach in the theory of strongly correlated electron systems. In 
such a case the methods of mean-field type are useful we have used one of these methods. The justification 

of an approximation is the reproduction of limiting cases and known exact results, and also the accordance with 
experimental data. Discussing the generalized mean-field approximation used by us we note the following facts. 

The approximation reproduces the exact atomic and band limits in the Hubbard model and also describes metal- 
insulator transition as mentioned above. For the Hubbard model the criterion of metal-insulator transition obtained 
by means of this approximation is Uc = 2wo {wq = z\to\) in agreement with the Mott's criterion [pT| and the 
result of the random dispersion approximation Q in infinite dimensions. In the limit of strong interactions the 
approach satisfies another requirement of a good approximate theory for the Hubbard model at half-filling formulated 
in Ref. = 2wo (W/, Wu are the widths of the lower and upper Hubbard bands). For the special case 



t' = of the model (2.2) the generalized mean- field approximation reproduces exact criterion (2.8) of metal-insulator 
transition, the ground state energy and doublon concentration. For very small (large) values of electron concentration 
and U ~^ oo the approximation gives the band description according to the general physical ideas. The approach 
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allows I^J] also to describe the observed transition in some Mott-Hubbard compounds from a metallic state to an 
insulating phase with the increase of temperature in a paramagnetic state (see, for example, book |]5l| , p. 178, 195). 

On the other hand, a deficiency of the generalized mean-field approximation is the breakdown of Fermi liquid 
behaviour of the systems with weak interactions and half-filled band (however, it should be noted that in the present 
paper we consider the cases of strong and intermediate interactions). The question about Fermi liquid behaviour 
(i.e. the fulfilment of Luttinger's theorem |^^) at intermediate to strong interactions in three dimensions is under 
discussion |^,^. This is caused by the reason that Luttinger's theorem is based [^,^ on the perturbation expansion, 
assuming an adiabatic relation between a noninteracting system and an interacting one, which may not be fulfiled 
for the case of intermediate and strong interactions. In particular, the authors have found the violation 

of Luttinger's theorem in the doped Hubbard model, and t-J model which is a partial case of the former one. 
Therefore, although the system is metallic away from half-filling, it may be non-Fermi liquid at intermediate and 
strong interactions because of a violation of Luttinger's theorem. 

Another limitation of the generalized mean-field theory is that for small values t' /U the approximaton does not 
reproduce the physics of the interaction of local magnetic moments as this holds in the framework of the effective 
Hamiltonian method. Nevertheless, we consider the paramagnetic phase, and in this case the effective inter-atomic 
exchange interaction (kinematic superexchange) leads to the renormalization of chemical potential only (in the molec- 
ular field approximation). Thus, apparently, it can be assumed that the generalized mean- field approximation used 
by us can provide a good description of effects of electron correlations in the paramagnetic phase and the region of 
intermediate (when a metal-insulator transition can occur) and strong interactions. Note also, that a modification ]5^ ] 
of this approximation proposed one of us allows to take into account effects caused by the e ffecti ve interaction of local 
magnetic moments; at the same time, the main features of quasiparticle energy spectrum (2.12) still persist. 

The noted generalized mean-field approximation has been applied also for study of metal-insulator transition and 
properties of generalized Hubbard model with correlated hopping [p5|^^8| , for description of metal- insulator transition 
in the doubly degenerate Hubbard model 12^,^ and a doubly orbitally degenerate model with correlated hopping . 



III. CHEMICAL POTENTIAL OF THE MODEL AT HALF-FILLING 



Consider the important situation of half-filled band n 
((XO) = (Xf)): 



1 . Chemical potential of the model is given by the equation 



E 



-Ei(k) , , ' E2(k) , -, 

exp + 1 exp + 1 



E 



exp 



--Ei{k) , T ' --E2(k) I 1 

TT^ + 1 exp TT—'- + 1 



(3.1) 



From Eq. (3.1) for the rectangular density of states and T = we find chemical potential of the generalized Hubbard 
model with correlated hopping in the region of metal-insulator transition: 



-U 



{U <w + w), 



w + w 

U w ~ w , ^ 
^=— H — (U>w + w), 



(3.2) 
(3.3) 



with w — z\t\, w — z \t\ beeing the half-widths of the lower and upper Hubbard band respectively. For the Hubbard 



model formulae (3.2) and (3.3) give the well-known result /i = U/2. At t' = chemical potential of the generalized 
Hubbard model is equal to n = U/2 beeing a consequence of the electron-hole symmetry which is a characteristic of 
the model in this case. Note that the value Uc = w -\- w corresponds to the metal-insulator transition point of the 
generalized Hubbard model with correlated hopping. 

At arbitrary value of temperature from Eq. (3.1) we find the expression for calculating chemical potential: 



1 



1 



exp 



-E2(e) 



exp 



de = 0, 



(3.4) 



where Ei{e), E2{e) are obtained from the respective formulae ( p.2C| ) for i?i(k) and -E2(k) substituting tk^e, ^k^fe, 
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Figs. |[^, where the dependences of chemical potential on the ratio U/w and temperature are plotted, show that 
chemical potential of the generalized Hubbard model with correlated hopping fi > [7/2 and only in the absen ce o f 
correlated hopping or within high temperature region the chemical potential is /it = t//2. From Eqs. (3.2) and (3.3) 
one can see that chemical potential of the generalized Hubbard model with correlated hopping becomes dependent on 
w, w. The dependences on these parameters are different in metallic and insulating phases, this leads to a kink 
at the point of metal-insulator transition (see Fig. |l|); < = (i.e. w = 0) and t' = O.bt correspond to the values of 
correlated hopping parameters ti — Ti/\to\ — and T2 — T2/\to\ = 0.5. Taking into consideration correlated hopping 
leads to the non-equivalence of the lower and upper Hubbard bands. This changes the value of chemical potential 
from the point fi = U/2: with the increase of the correlated hopping T2 (i.e. with decreasing w) chemical potential 
moves towards the upper Hubbard band. 

With the increase of temperature the kink in the chemical potential curve disappears (Fig. ^). Similar peculiarity 
in the evolution of free energy dependence on parameter U/w with temperature was noted by Mott [5l| ]. 

From Fig. |^ one can see that in the model under consideration within low and room temperature regions chemical 
potential is essentially dependent not only on the parameters w and w, but becomes also dependent on temperature (in 
contrast to chemical potential of the Hubbard model). Moreover, with the decrease of temperature chemical potential 
rapidly increases, this dependence is different for the different values of correlated hopping parameters ti, T2. The 
found temperature dependence of chemical potential can be explained by the following reasons. At zero temperature 
in an insulating phase the chemical potential is located at the centre of the energy gap (between the top of the lower 
Hubbard band and the bottom of the upper Hubbard band), the chemical potential value is larger than U/2. As the 
temperature increases some electrons are thermally activated and the energy levels within the upper (lower) Hubbard 
band become to be occupied (empty). The width of the lower Hubbard band is larger than the width of the upper 
Hubbard band, and thus, the density of states in the lower band is smaller than the density of states in the upper band 
(as a result that the distance between the nearest energy levels within the lower Hubbard band is larger versus the 
corresponding distance in the upper Hubbard band). Consequently, with increasing temperature chemical potential 
position moves towards the lower Hubbard band, and chemical potential value decreases. 

Within high temperature region in the generalized Hubbard model chemical potential tends to U/2 with the increase 
of temperature; really, for T — » 00 the probabilities of an electron finding within the lower and upper Hubbard bands 
are equal. 



IV. CHEMICAL POTENTIAL OF THE DOPED MOTT-HUBBARD SYSTEMS 



Find the dependence of chemical potential /i of the generalized Hubbard model with correlated hopping on carrier 
concentration. Chemical potential is given by the equation {{Xf)-{Xf)=n-l): 



n 

2N 



1 ^n/2 
N 



exp 

E 



1 exp 
^k 



exp 



^ + 1 exp^-t-1 



n — 1, 



with Ei(k) and i?2(k) determining by expression (2.12), 



g(k) - e(k) + U 
2Q(k) 



1 g(k)-£(k) + [/ 

2 2g(k) 



Q(k) is def ined by formula ( p.l3| ), and e(k), e(k) are given by expressions ( 2.14 ). 
Eq. (4.1) at T = is rewritten in the form: 



N 



^^[i?k(n-l) + l-n/2]0(i?2(k)) 

k 

I- J2 [^k(l -n) + 1/2] eiE^ik)) = 1 - n/2. 



(4.1) 



(4.2) 



(4.3) 



We assume the rectangular density of states and pass in Eq. (4.3) from a summation over k to an integration over 
energy. We solve this equation for the case of a paramagnetic phase of the doped Mott-Hubbard insulator (the energy 
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spectrum is supposed to has form ( ^.1^ ), ( 2.19| )), and find that chemical potential of the generalized Hubbard model 



(3n-2)(2-2n + n2) 

M = ^ {2-nY ^'"^'^^ ^" ^ ^' ^ ^ ^ 

(3n-4)(2-2n + 7i2) 

Ai = C/+^ '-^^ -w{n) (n>l), (4.5) 



with w(n) = wo(l — nri), ?D(n) = wo(l — nri — 2t2). 



Concentration dependence of the chemical potential (4^) ancH4.5) is plotted in Fig. For the Hubbard model this 



dependence is in qualitative agreement with the results of Refs. [p9| , |60[ , obtained by means of the Gutzwiller variational 
method |l6l|] in three dimensions and the projective self-consistent technique in infinite dimensions. However, the 
curve camber in Fig. || is as opposed to that plotted in paper the second derivatives of chemical potential 

with electron concentration n are opposite in sign; thus, the behaviour of electron compressibility of system (which is 



defined by the quantity ^), determining with the help of Eq. (4.4) and (|4.5| ), is inverse to the compressibility behaviour 



obtained using /i(n)-dependence of Rcf. [Q. The results of works [^|60| coincide with the conclusions derived from 
Fig. U, however, there is another essential difference between them (see the next paragraph). It should be also noted 
that our calculations show ^ > in all region of the electron concentration. Thus we do not predict the existance 
of phase separation upon doping in the framework of the model under consideration in agreement with the quantum 



Monte-Carlo results |62| for the Hubbard model in infinite dimensions. Such an instability is experimentally observed in 
high-temperature superconductors and manganites with colossal magnetoresistance [ p4[ . This interesting feature 
was predicted theoretically to be a most common behaviour of extended Hubbard model with nearest-neighbour 
Coulomb repulsion and double exchange model [ |65| , t-J model |63| , and pseudospin-electron model JSGf . 

Fig. ^ displays the jump in chemical potential at half-filling in going from electron to hole doping, this discontinuity 
in chemical potential is strictly equal to the energy gap width /S.E = U — w — w. Therefore, our calculations confirm the 
first scenario of the chemical potential evolution upon doping mentioned in Section 1. As a result of the absence of a 
rigorous result on the chemical potential behaviour, one of the main tests of an approximate theory is the comparison 
with experimental data. One knows no experimental data, which show that the theoretical prediction of some authors, 
such as the existence of mid-gap states in the Mott-Hubbard insulator upon doping, indeed occur in real transition 
metal compounds. However, the experiment ||6^ , |68|| suggests the jump in Fermi energy near half-filling is equal to the 
Mott-Hubbard gap in some strongly correlated electron systems. 

Fig. ^ shows that correlated hopping changes essentially the position of chemical potential of the generalized 
Hubbard model: at n < 4/3 increasing correlated hopping Ti (or T2) leads to the increase of chemical potential, and 
for 71 > 4/3 chemical potential decreases with increasing correlated hoppings Ti, T2. This can be explained by the 
following arguments. For the electron concentration n — > 1+ the chemical potential is located at the bottom of the 
upper Hubbard band. When the electron concentration is n ^ 2^, the chemical potential position corresponds to 
the top of the upper Hubbard band. For the electron concentration n — > 1"*" chemical potential of the generalized 
Hubbard model is larger than chemical potential of the Hubbard model, and vise versa for the electron concentration 
n — > 2~, as a result of the fact that width of the upper band of the Hubbard model is larger than width of the upper 
band of the model with electron- hole asymmetry. Thus, the value uq of electron concentration must exist satisfying 
the requirements: at this value the chemical potentials of the Hubbard and generalized Hubbard models are equal 
(in this case the chemical potential is located at the centre of the upper Hubbard band, i.e. fi ~ U), for electron 
concentration n < uq increasing correlated hoppings Ti, T2 leads to the increase of chemical potential, and for n > uq 
chemical potential decreases with increasing correlated hoppings Ti , T2. 

Similarly, we can explain the dependence of chemical potential of the generalized Hubbard model on the correlated 
hopping Ti for electron concentration n < 1. But in this case, with contrast to the previous situation, the dependence of 
chemical potential on the correlated hopping Ti is weak (besides, for the doped Mott-Hubbard insulators with n < 1 
chemical potential has no dependence on the correlated hopping T2). It should be stressed that the found values 
no = 2/3, 4/3 are the electron concentrations corresponding to the change of conductivity type of the system p^ . 

From Figs. ^, 0one can see that for electron concentrations n > 1 the chemical potential of the generalized Hubbard 
model has the linear dependence on correlated hopping parameters, and this dependence on T2 is stronger versus ti. 
For the values of correlated hopping parameters ri =0, T2 = 0.5 the chemical potential of the generalized Hubbard 
model at n > 1 becomes equal to the chemical potential of the Hubbard model in the atomic limit fi = U (though 
to 7^ 0); in this case chemical potential of the generalized Hubbard model does not depend on electron concentration. 
Notice that for these or close values of correlated hopping parameters within the model under consideration the 
following situation, which is analogous one observed in heavy fermion substances and compounds with mixed valency, 
can occur. Near the Fermi level which lies within the lower Hubbard band, the level or very narrow upper Hubbard 
band (being analogous to the /-level in heavy fermion substances and compounds with mixed valency) can arise in 
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conduction electron subsystem, as a result of electron-electron interactions of correlated hopping type. In particular, 
this can lead to the very strong enhancement of effective mass of carriers 27 2^. 

It should be noted that using more realistic expressions for the density of states does not change essentially the 
behaviour character of the obtained dependences of chemical potential on the model parameters. For example, in the 
case of semielliptic density of states: 



2 



we obtain the following equation for chemical potential of the generalized Hubbard model with correlated hopping: 

(4.7) 




3n 1 

T " 2 



{2-2n + n'^)w{n) 
{fi - C/)n2 r 7 (/i -U)n y 3n 



27r(2 - 2n + n2)w(n) y \{2 - 2n + n'^)w{n) J 4 



(n< 1), 

(4.8) 

{n > 1). 



Analysis of these equations shows that the obtained behaviour peculiarities of chemical potential for the case of the 
rectangular density of states persist. Moreover, a choice of the semielliptic density of states does not change even the 
values no = 2/3, 4/3. 



A. Temperature dependence of chemical potential 

Find the dependence of chemical potential of the generalized Hubb ard model with correlated hopping upon tem- 
perature. Using the rectangular density of states, we rewrite Eq. ([4.1|) for chemical potential in the form: 

- 1. (4.9) 



(n) w{n) 

de 1 — n/2 f de 



J exp^-fl 2?«(n) J exp + 1 

After integration for the case of the doped Mott-Hubbard insulator we obtain the following equation for chemical 
potential: 

1 1 + CXD „ /2 _ 1 1 + CXD 

In ^ P ' , ^gg-iln ' + "^P , (n<l), (4.10) 



\w{n) 1 + exp 2iiw{n) l + cxp 



1 1 + exp -A'+C/-/3A(n) ^ 1 exp A'+(2-»)^(n) 



In hrrr = In , I , {n > 1), (4.11) 



4/3w(n) 1 + exp -t^+ufMn) 4w{n) 1 + exp 



with the coefficients (3 = ^~^:^+" ^ - 2-2n+n 



JO. 



On the base of expressions ( I.IC ) and ( 4.11 ) we obtain the dependence of chemical potential on electron concentra- 
tion at some values of temperature (Fig. |), and //(T)-dependence for some values of electron concentration (Fig. ||). 
Here we limit ourselves to a consideration of the Hubbard model case, because taking into account correlated hopping 
of electrons at T 7^ does not lead to qualitatively new distinctions of the chemical potential behaviour from the 
mentioned one at T = 0. 

From Fig. ^ one can see that at temperature T > the discontinuous jump in chemical potential, obtained in the 
previous subsection, becomes more smooth, the dependence of chemical potential on electron concentration becomes 
a S-shaped curve with flex point at the electron concentration n — I. 

Figs. ^, ^ show that for the electron concentration n < no chemical potential of the doped Mott-Hubbard insulator 
decreases with the increase of temperature, and for n > hq the increase of temperature leads to the increase of 
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chemical potential; at n = no the chemical potential does not depend on temperature. These peculiarities of the 
temperature-dependence of chemical potential can be explained by the following reasons. Let us consider the case 
n < 1. For small values of electron concentration n < Uq the Fermi energy lies below the centre of the lower Hubbard 
band, thus with the increase of temperature electrons become to occupy the above-lied empty energy levels, this leads 
to decreasing chemical potential (an analogue of the case of free electrons). For the electron concentration n — hq 
the Fermi energy is located at the centre of the lower Hubbard band, as the temperature increases the above-lied 
empty energy levels within the lower Hubbard band become to be occupied, but the chemical potential does not 
change its value as a result of symmetry about the centre of the band. For enough high temperatures and values 
of the intra-atomic Coulomb repulsion U > 2wo some electrons are thermally activated into the upper Hubbard 
band, therefore the chemical potential value decreases with increasing temperature at n = hq. For values of electron 
concentration n > hq the number of empty energy levels within the lower Hubbard band is small because the Fermi 
energy in this case lies above the centre of the lower band, consequently with the increase of temperature chemical 
potential increases owing to occupancy of these levels and a presence of the empty upper Hubbard band. Similarly, 
the influence of temperature and electron concentration on chemical potential of the Hubbard model for the values of 
electron concentration n > 1 can be explained. 

Note that one can see an analogy of the temperature dependence of chemical potential of the model under con- 
sideration with the same dependence in doped semiconductors. Really, for electron concentration n < uq within the 
model the system has n-type of conductivity and chemical potential decreases with increasing temperature as 
in 71-type semiconductors. For electron concentration n > hq the system has p-type of conductivity, correspondingly 
chemical potential increases with the increase of temperature, as it is characteristically for p-type semiconductors. 

V. CONCLUSIONS 

In the present paper chemical potential of the generalized Hubbard model with correlated hopping has been studied 
by means of a generalized mean-field approximation. The peculiarity of the model in comparison with similar gener- 
alized Hubbard models is the concentration dependence of hopping integrals caused by taking into consideration the 
matrix elements of electron-electron interaction which describe inter-site hopping of electrons (correlated hopping) . 

The chemical potential as a function of the model energy parameters (U, t{n), t(n), f{n)), electron concentration 
and temperature has been found. At half-filling, in contrast with the Hubbard model where /i = U/2, the dependences 
of chemical potential on energy parameters of the model are different in metallic and insulating phases leading to a 
kink at the point of metal-insulator transition at zero temperature; with the increase of temperature the kink in the 
chemical potential curve disappears. With the increase of the correlated hopping T2 (i.e. with decreasing w) chemical 
potential moves towards the upper Hubbard band. 

Chemical potential of the generalized Hubbard model with correlated hoping is temperature-dependent even in the 
case of half- filled band (in contrast to the Hubbard model). With the decrease of temperature chemical potential 
rapidly increases, this dependence is different for the different values of correlated hopping parameters ri, T2. In high 
temperature region chemical potential of the asymmetric Hubbard model tends to U/2 as temperature increases. 

In strong coupling regime we have found that correlated hopping changes essentially the location of chemical 
potential of the generalized Hubbard model: for the electron concentration n < uq increasing correlated hopping 
parameters leads to the increase of chemical potential, and for the case n > uq chemical potential decreases with 
increasing correlated hopping parameters. Similarly, for the electron concentration n < uq chemical potential decreases 
with the increase of temperature, and for n > hq the increase of temperature leads to the increase of chemical 
potential; at the value of electron concentration n = uq the chemical potential does not depend on correlated hopping 
and temperature. The found values no — 2/3, 4/3 correspond to the electron concentration when conductivity type 
of the system changes; chemical potential of the systems with these electron concentrations is located at the centre 
of the lower and upper Hubbard band respectively. 

At zero temperature we have found the jump in chemical potential at half-filling in going from electron to hole 
doping, this discontinuity in chemical potential is strictly equal to the energy gap width AE = U — w — w. If 
the temperature is greater than zero then the discontinuous jump of chemical potential becomes more smooth, the 
dependence of chemical potential on electron concentration becomes a S-shaped curve with flex point at the electron 
concentration n = 1. 
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FIG. 1. Dependence of chemical potential ji of the 
generalized Hubbard model with correlated hopping 
in the ground state: the upper curve corresponds to 
Ti = 0, T2 = 0.5; the middle curve - ti = 0, T2 = 0.3; 
the lower curve - n = T2 = (the Hubbard model). 
The asterisks denote the point of metal-insulator 
transition. 
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FIG. 2. Evolution of the dependence of chemi- 
cal potential /x of the generalized Hubbard model 
(ti — 0, T2 = 0.5) on U/w with temperature 
e/w = kBT/w = 0, 0.1, 0.3, 1 respectively. The low- 
ermost curve corresponds to the behaviour of chemi- 
cal potential of the Hubbard model. 




FIG. 3. Temperature dependence of chemical po- FIG. 4. Chemical potential of the generalized Hub- 

tcntial /i of the generalized Hubbard model with bard model versus the electron concentration n at 

correlated hopping at U /2w = 1: the upper curve U/wq =4: n = T2 = (1); n = T2 = 0.2 (2); 

corresponds to n = T2 = 0.3; the lower curve - n = 0, r2 = 0.5 (3). 



Ti = T2 = 0.1; the straight - Ti = r2 = (the Hub- 
bard model). 



12 



6.00 



4.00 



~2.00 



-0.00 



-2.00 




0.00 



2.00 



FIG. 5. Evolution of the concentration de- 
pendence of chemical potential of the Hub- 
bard model with temperature for U/wo = 4: 
9/wo = knT/wo = 0, 0.1, 0.3, 0.5 respectively. 
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FIG. 6. Chemical potential of the generalized Hub- 
bard model as a function of the correlated hopping pa- 
rameter n = Ti/ltol for U/wo = 4: the lower curve 
corresponds to n = 1.7, the middle - n = 4/3, the 
upper - n = 1.2. 
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FIG. 7. Chemical potential of the generalized Hub- 
bard model as a function of the correlated hopping pa- 
rameter T2 = T2/\to\ for U/wo = 4: the lower curve 
corresponds to n = 1.7, the middle - n = 4/3, the 
upper - n = 1.2. 
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FIG. 8. Temperature dependence of the chemical 
potential of the Hubbard model with U/wo = 8 
for different values of electron concentration: 
n = 0.98, 0.9, 2/3, 0.5, 0.3 (from top to bottom). 
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